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Abstract 

In this paper, we present a model describing the time evolution of two dimensional surface waves 
in gravity and infinite depth. The model of six interacting modes derives from the normal form 
of the system describing the dynamics of surface waves and is governed by a Hamiltonian system 
of equations of cubic order in the amplitudes of the waves. We derive a Hamiltonian system with 
two degrees of freedom from this Hamiltonian using conserved quantities. The interactions are 
those of two coupled Benjamin-Feir resonances. The temporal evolution of the amplitude of the 
different modes is described according to the parameters of the system. In particular, we study 
the energy exchange produced by the modulations of the amplitudes of the modes. The evolution 
of the modes reveals a chaotic dynamics. 
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I. INTRODUCTION 

In 1967, Benjamin and Feir studied the linear stability of periodic solutions of per- 
manent form for gravity surface waves in infinite depth. They showed that wave trains are 
unstable to modulation perturbation i.e that instabilities occur due to the interaction of the 
carrier wave and two wave perturbations. This kind of instability has been first found by 
Lighthill [2] and confirmed theoretically and experimentally by Benjamin and Feir. This 
result was also confirmed by Zakharov [3] using a Hamiltonian formalism for the surface 
wave problem. This instability was interpreted in terms of resonance conditions between 
the carrier wave and the two side bands. Later in 1977, Lake et al. j^J observed in wave 
tank experiments the first modulation and demodulation of the recurrence Fermi-Pasta- 
Ulam phenomenon 

Many works using the Hamiltonian formalism which was first derived by Zakharov, have 
been carried out to understand the time evolution of the resonant interaction phenomenon 
Q, 0| . We refer in particular, for three dimensions, to Shemer and Stiassnie p , Stiassnie 



and Shemer llCI, Shrira et al. Illfl and Badulin et al. 



121. In two dimensions, we refer to 



[nlshrira et al. Q 

Benzekri et al. 13], Caponi et al. [14f who use a seven modes model in order to show a 

n 

chaotic behavior. In order to study the long time evolution of surface waves, Zufiria [15| 
derived from the Zakharov water waves equations a simple model of three interacting modes 
[Eq.(lO) in Ref. [15]] which leads to a chaotic behavior. However this model is not invariant 
by time reversing symmetry. 

Other works using the Hamiltonian formalism in fluid interfaces were carried out in 
Refs. 

In order to give a qualitative description of the spatio-temporal evolution of the dynamical 
system described by a Hamiltonian water wave system, we propose in our work a simple 
model using all the symmetries of the initial equations of surface waves. In particular, we 
use the time reversing symmetry to derive a system of Hamilton's equations describing the 
dynamics of Benjamin- Feir resonant waves. This symmetry plays an important role in the 
study of the system. The evolution of the complex amplitudes of the modes is governed by 
a cubic dynamical system of equations. The system is then a set of six ordinary differential 
equations. We show the existence of invariants that reduces this system to two degrees of 
freedom. 
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We use a geometrical approach to understand the energy exchange mechanism among the 
modes of the wave. For some initial conditions, we find that the Benjamin-Feir resonance 
leads to a continual periodic energy transfer from one component of the wave to the others. 
We also find by numerical simulations an irregular transfer of energy among the modes of 
the wave. The energy transfer of this kind is related to the existence of chaos in our system. 

This paper is organized as follows: In Sec. |Hj we describe the model which is an approx- 
imation of the Zakharov Hamiltonian water waves system developed by Krazitskii [18]. In 
this model, the interactions require the coexistence of two resonating four-wave interactions. 
The Hamiltonian of this model written in normal form is of order four in the amplitudes of 
the wave. 

The system has six degrees of freedom, then having six equations for the complex ampli- 
tudes. We find four integrals of motion thanks to the Manley-Rowe relations. Therefore, we 
reduce the model to a Hamiltonian system with two degrees of freedom. 

The knowledge of these invariants allows us to subdivide the phase space in different 
regions and determine analytical solutions like traveling and standing waves. In particular, 
we analyze how the system behaves starting from different initial conditions. Section IIHI is 
devoted to the study of the dynamics of the system with only one resonance. The study of 
the chaotic dynamics with the two interacting resonances is done in Sec. IIVI 



II. THE MODEL 



In an inviscid and incompressible fluid approximation, the equations of the two- 
dimensional gravity waves on 

A<f(x,y,t) — 0, for — oo < y < rj(x, t), (1) 

gr ] (x,t) + ^(x,y,t) + ^\\7ip(x,y,t)\ 2 = 0, for y = rj(x, t), (2) 
dt] dip dt] dip 

— (x, y, t) + — (x, y, t)-Q^{x } V, t) - -q-(x, V, t) = 0, for y = T](x, t), (3) 

lim Vip(x,y,t) = 0, (4) 

y— oo 

where the unknown variables are the surface shape of the wave rj, and the stream function 
ip and g is the strength of the gravity. Zakharov [3J showed that these equations can be 
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expressed in a Hamiltonian form: 

da k OH ^ 

^ = -^, fa^M, (5) 

where the complex amplitudes a k are linear combinations of the Fourier coefficients of i] (x, t) 
and of the free surface velocity potential ifi(x,t) defined by 

ip(x,t) = <f(x,T](x,t),t). 

If we restrict ourself to the spatially periodic waves then fcel More precisely, if we define 
the Fourier expansion of the wave 

^ +00 



T](x,t) = — J2 Vk 

III , 

k=—oo 



Akx 



the relation between the coefficients r\ k and the complex symplectic coordinates a k is 

Vk(t) = —^-\k\ 1/4 (a k + a*_ k ) . 

The Hamiltonian is given by the total energy of the system when expressed as function of 
the a k (see Ref. 

Our model is obtained by truncating the general (infinite dimensional) Hamiltonian de- 
rived by Zakharov in the form stated by Krazitszkii [18] up to order four and retaining a 
minimal set of modes in order to describe a Benjamin- Feir instability for gravity surface 
waves in infinite depth. 

In order to preserve time reversal symmetry, this minimal number of complex modes is six, 
±ki, ±k 2 , ±/c 3 that leads to a set of two coupled resonances verifying: 

kx + k 2 = 2k 3 , and u kl + u k2 = 2u k3 , (6) 

together with the dispersion relation: 



Uk = yg\k\, for k = ±kx,±k 2 ,±k 3 . (7) 

Dyachencho and Zakharov 19| described the general families of such type of resonances. 
In our case this leads to only one family of resonances, namely 

ki = —a, k 2 = 9a, k 3 = 4a, for aeZ*. 
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For the simplicity of the exposition, the modes will be denoted by aj instead of for 
j = 1, 2, 3 and a_j instead of a_&.. 

It is clear from Eq. (JHJ) that such resonances are built from central carriers, a 3 and a_3, and 
two equally distant side bands. 

The expression of the Hamiltonian in the variables aj up to order 4 is given in the Appendix. 
Then following Ref . [18[ , we perform a canonical transformation from the variables aj to the 
variables bj such that the Hamiltonian is mapped into its normal form. We rederive it in 
the Appendix. 

The Hamiltonian of the system becomes in the variables bj 

n= E ^-(N 2 + iM 2 )+ E ^(N 4 + M 4 -4N 2 M 2 ) 

3=1,2,3 4=1,2,3 

+v 21 (M 2 N 2 + \b-2\ 2 \b-i\ 2 ) + v_ 21 (N>-i| 2 + M 2 N 2 ) 
+v 3l (-N 2 N 2 - M 2 M 2 ) + K-31 (N 2 M 2 + IM 2 N 2 ) 

+V 23 (\b 2 \ 2 \b 3 \ 2 + |6_ 2 | 2 |&_ 3 | 2 ) + ^-23 (-|& 2 | 2 |&-3| 2 - I&-2H&3I 2 ) 
+iiy(6!6 2 63 2 + b\b* 2 b\ +b_ 1 b_ 2 b* 2 3 + b*_ 1 b*_ 2 b 2 _ 3 ), (8) 

where the coefficients Uj, and W are positive. The values of all these coefficients are 

given in the Appendix . 

The equations of motion are : 

dbj_ .dH 
~dt ~ ^dbj' 

for j = -1,-2, -3, 1,2, 3. 

We notice that, in the infinite dimensional showed by Dyachenko and Za- 

kharov [19] and Craig and Worfolk j^J, the canonical transformation that eliminates 
the cubic terms also cancels the fourth order terms. In our case, starting from the finite 
dimensional truncation, the Benjamin- Feir fourth order terms are instead present in the 
normal form. This shows a difference between the infinite dimensional case and our finite 
dimensional model. A similar Hamiltonian model for fluid interfaces has been derived in 
Ref. where Benjamin- Feir resonances appear even in the infinite dimensional case due 
to the presence of the non-vanishing density of the upper fluid. 

Using some scaling relations, it is possible to assume that u>i = 1, u 2 = 3 and u 3 = 2. 
Indeed, if we rescale the variables bj by a factor \f\ for A 6 M* + , i.e., if we replace the 
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Hamiltonian H(bj,b*) by XH (bj / A/A, b* j \/A) which is a transformation that preserves the 
equations of motion, and if we rescale time by the same factor A, i.e., we multiply the 
Hamiltonian by A, then it is easy to see that the quartic part of the Hamiltonian is un- 
changed and the quadratic part is multiplied by A. By choosing A = 1/ \Jg\a\, it means that 
we can consider that the frequencies are lo\ — 1, uj 2 — 3 and oj 3 = 2 without loss of generality. 

As we shall see, the existence of several integrals of the motion in our system, despite 
the fact that they are not sufficient to turn it into an integrable one, allows to reduce the 
number of degrees of freedom from six to two and to have important information about the 
dynamics in phase space. 

The equations of evolution of |&i| 2 , \b 2 \ 2 and |6 3 | 2 are 



> 2 


= «y( 6 M - 




d \h I 2 


= t—(bib 2 b 3 - 


blb*bl), 


-N 2 


= -iWihbzbf 





Therefore we see that the following two quantities are conserved 

1+ = 2|6 1 | 2 + |6 3 r, (9) 
J+ = |6i| 2 + 3|6 2 | 2 + 2|6 3 | 2 . (10) 

In what follows, we will also use a linear combination of these two conserved quantities 

K + = N 2 -N 2 = i(2I + -J + ). (11) 

By defining Nj = \bj\ 2 for j = 1, 2, 3, the constants of motion become: 

1+ = 2N, + N 3 , 

J+ = JVx + 3N 2 + 2N 3 . 

Similarly, we have two additional conserved quantities for the other resonance 

I_ = 2N_! + iV_ 3 , 

J- = iV_i + 3AL 2 + 2iV_ 3 , 

where N_j = \b^\ 2 for j = 1, 2, 3. 

It can be checked that the four integrals I + , J + , I_, J_ are independent and in involution. 
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The key point to be noticed here, is the fact that each one of these first integrals depends 
on the variables of one resonance alone even if, as it will be seen, there is a dynamical 
interaction of the two sets of variables. 

Let us now look at the structure of the invariant surfaces in phase space in view of the 
previous integrals of motion. Even if they can not completely fix the orbits of the system, 
they help to understand many properties of the different dynamical regimes observed in 
this model. 

In the space E + = (Ni, N 2 , N%) (respectively E_ = (iV_i, AL 2 , iV_ 3 )), J + = const, (resp. 
J_ = const.) defines a portion of a plane, that we represent in Fig. ^ On the other hand, 
the intersection of this set with I + = const, (resp. I_ = const.) defines a family of segments 
C + (I + ) (resp. C_(I_)). 

Let us consider the structure of E + and let us fix J + 7^ 0. For I + = 0, the set C + (0) 
reduces to one point P whose coordinates are P = (0, J+/3, 0). 

For < 1+ < the intersection of the two invariant surfaces in E + , C+(I+) is a segment 
close to P. 

When I + = J+/2, C+(JJ+/2) is a segment with one extremity T = (0,0, J+/2). We notice 
that on this curve we have N\ = N 2 . 

For J + /2 < I + < 2J + , the set C+(I+) is a segment close to S = (J + , 0, 0). 

For I + = 2J + , the intersection is the point S. 

For I + > 2J + or for I + < the intersection is empty. 

We notice that E + x E_ needs to be supplemented by additional variables like the six 
angles $j of the complex variables bj in order to build the phase space. Therefore, since 
the projection of an orbit on E + x E_ alone is only part of the information we must look 
at it as a partial, although useful, information on the dynamics. This is why, in each case 
we shall supplement its description with additional information in order to make a complete 
characterization of the dynamics. 
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III. DYNAMICS WITH ONE RESONANCE 

To get insight into the full dynamics, we begin with simple and integrable situations 
which will allow us to understand the general case with two resonances. We start exploring 
the phase space by choosing initial conditions leading to simple situations for which we can 
exhibit different kinds of shapes of the wave. A complete description of phase space is given 
by a numerical analysis. 

It follows from the conserved quantities, that an orbit corresponding to an initial condition 
with all null coordinates in one of the two resonance variables Nj will keep this property for 
any time. For instance, we assume that J_ = 0. This implies that iV_i = N_ 2 — N_ 3 = for 
all time. Thus, we first study the dynamics with only one resonance which is an integrable 
case since Hamiltonian (fT2"j) has three degrees of freedom and three conserved quantities (e.g. 
H, I + , K + ). Hamiltonian JHJ) describing the dynamics in the variables (bi,b 2 ,b 3 ) becomes : 

H = fi^ 2 + 3\b 2 \ 2 + 2\h\ 2 + UM 4 + U 2 \b 2 \ 4 + U 3 \b 3 \ A 

+^ 2 iN 2 |6i| 2 - v 31 \b 3 \ 2 \h\ 2 + v 23 \b 2 \ 2 \b 3 \ 2 

+\w(b 1 b 2 bf + b* 1 b* 2 b 2 3 ). (12) 

In what follows, we describe the dynamics on the sets C+(I+) described in the previous 
section (see Fig. [I}. We start with the simplest cases which are the three points P, T and 5* 
and the segment {PS) (one of the boundaries of the triangle E + ) for which we can exhibit 
analytical expressions for the shape of the wave i](x,t). 

A. Dynamics of orbits of P,T,S 

First we consider the orbits corresponding to the point P of E + (where I + = 0). The 
equations of motion read 



iVi EE 


N 2 


= o, 


N 3 = 


N 2 


= o, 


N 2 EE 


N 2 


_ J+ 



s 



^| = (3 + 2f/ 2 N 2 )& 2 . 
The solution in the variables bj is then 

b 2 (t) = ^exp(-i(3 + 2[/ 2 ^±)t-i$ 
bj(t) = 0, for j = 1,3. (13) 

Therefore for any J + > these orbits are periodic in time. Since we have only one resonance, 
the original variables aj are equal to bj (see Remark 1 of the Appendix). 
Thus the expression of the shape of the wave is: 

fo —1/4 

V (x, t) = V 9 (b 2 (t) exp(9ix) + b* 2 (t) exp(-9zz)) , (14) 

Z7TV 2 

which is also 

77(37, t) = 77 cos(9(x - c 2 t) + 0), (15) 
i.e. traveling waves with velocity 

1 471-2 

c2 = - + —V9U 2 vo ■ (is) 

A plot of the shape of this wave is represented on Fig. El 
Similarly for S, we obtain traveling waves running backwards t]{x,t) = ?7ocos(x + c\t + 0) 
with velocity 

Cl = l + 4<K 2 ^jU 1 r ]0 2 . 

For T, we obtain traveling waves running forwards rj(x,t) = t] Q cos(4(x — c^t) + 0) with 
velocity 

c 3 = 2 + 71-2 'V9U3V0 2 ■ 

A linear stability of the points P,5* and T shows that these fixed points of the dynamics 
are elliptic. It means that nearby trajectories will mimic the behavior of the dynamics 
of these points. In particular, the surface waves near these points are time quasiperiodic 
waves close to traveling waves. 
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B. Dynamics on a point on the segment (PS) 

We consider the point which is at the intersection of C+(I+) and the segment (PS), for 
a value of I + G [0, 2JJ + ]. This point is such that 63(0) = and dbs(0)/dt = 0. From the 
equation for & 3 (t) , we deduce that bs(t) = for any time. From the two invariants I + and JJ + , 
we have that N\ and N 2 are constant and equal to Ni(t) = I+/2 and N 2 (t) = (2J+ — 1+)/6. 
The equations of motion become: 



where c x = 1 + (U x - V 21 /6)I + + V 21 /3J+ and c 2 = [3 + {-U 2 /3 + V 21 )l+ + 2/3V 2 J + ] /9. 
By taking into account Remark 1 of the Appendix, we have <%, = bj for j = 1, 2, 3. The 
expression of the shape of the wave is: 



which is the sum of two traveling waves with two different velocities. 

We notice that on P (respectively on S), the amplitude of the traveling wave with velocity 
Ci (respectively c 2 ) vanishes. Therefore the wave associated with P (respectively S) 
becomes the traveling wave with velocity c 2 (respectively C\) obtained in Sec. IIHI 




We deduce the expression of 61 and b 2 : 
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C. Dynamics of orbits on C+(I+) 



Next, we extend the dynamics with one resonance, by studying the dynamics on the 
full sets C + (I + ) for I + G [0,2J + ] . We determine numerically the profiles of the waves 
which are similar to those obtained previously namely stationary, traveling, time periodic 
and quasiperiodic waves. This study gives also all the information about the exchanges of 
energy which occur between the modes. 

In order to study the orbits in this case, since we have a Hamiltonian with three degrees of 
freedom and three constants of motion, we want to perform a canonical change of variables 
reducing the dimension of the system to one degree of freedom. Furthermore, we would 
like to change canonically the coordinates into action- angle variables (Nj,3>j) such that 



First we notice that if N 3 (0) = 0, then N 3 (t) = at any time t (see Sec. IIII B|) . In what 
follows, we always exclude this point which is at the intersection of {PS) and C + (I + ). 
In the case where I + = JJ+/2, we notice that since K + = Aq — N 2 is a conserved quantity 
(equal to zero in this case), if Aq vanishes at some time, then N 2 also vanishes. This 
dynamics occurs at the point T (see Sec. 1111 A|) . Except this particular case, there are two 
interesting cases where I + > J+/2 and I + < J+/2. If I + < J+/2, only Aq vanishes since 
N 2 > (J+ - 2I + )/3 > 0. If 1+ > J+/2, only iV 2 vanishes since Aq > (21+ - 2J + )/3 > 0. 

We remark that Aq is such that max(0,K + ) < Ni < I+/2, so the accessible phase space 
in the (pi, q\) = (Re(bi), Im(bi))-plane is a circle for K + < and an annulus for K + > 0. 
In order to avoid the two cases Aq = or A2 = for which the angle <p\ or 02 is not 
defined in the change of variables, we first show that there is only one trajectory Ti such 
that Aq(t) = 0, i.e which passes by zero in the (pi, q{)- plane. 

First we notice that Hamilton's equations associated with Hamiltonian (|T2*|) are invariant by 
the two transformations: 
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and 



b\e 



—it 



b 




where <\> is constant. Using these two invariants, one can show that there is basically one 
trajectory T\ which passes by zero in the (pi, gi)-plane. The closed trajectory T\ divides the 
phase space in the (pi, gi)-plane into two parts: the inside and the outside of T\. 
Similarly, for K + > 0, there is only one trajectory T 2 which passes by zero in (p 2 , q 2 ) -plane. 
We notice that N 2 = corresponds to N\ = K + and so a trajectory which passes by zero in 
the {p 2 , q 2 )-plenae, is a trajectory which passes by the point Ni = K + in the (pi, gi) -plane, 
which is the inside border of the accessible region. 

The trajectory T 2 divides the annulus into two parts: the inside and the outside of T 2 . Thus 
in any of these two parts, iV\ and N 2 are strictly positive. 

In the regions where all actions N\, N 2 and ^3 are strictly positive, we perform a canonical 
transformation that expresses Hamiltonian (fT2|) into action-angle variables : 



The actions N±, N 2 , N3 are respectively conjugate to the angles $1, $2, $3- Hamiltonian 
(JI2J) becomes 




(17) 



This transformation is canonical i.e. : 



idb j Adb* = > dN j A 




H 



iVi + 3N 2 + 2N 3 + U x Nl + U 2 Nl + U 3 N% 



+V 21 N 1 N 2 - Vs^Nz + V 23 N 2 N 3 




(18) 



Hamilton's equations for Hamiltonian (j!8)l are given by 



_ OH 



dt dNj ' 
dNj dH 



dt ' 9$,' 
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for j = 1,2, 3. 



Using the invariants I + and J + , we now proceed to the reduction of the dimension of the 
system to one degree of freedom. 

We perform the following canonical transformation from (iV\, N 2 , N 3 , $ 1; $ 2 , $3) to 
(N l ,K + ,I+,e 1 ,9 2 ,9 3 ): 

Ni = N h 
K + = Aq - N 2 , 
1+ = 2Aq + N 3 , 
61 = $1 + $ 2 - 2$ 3 , 

9 2 = -$2, 

e 3 = $ 3 . (19) 

The actions Aq, K + and I + are now respectively conjugate to angles 6\, 62 and 9 3 . In what 
follows, we use Aq instead of Aq for simplicity. Hamiltonian (fTSJl becomes 

H = n + QiNi + Vt 2 Nl + W^N~ 1 ^N 1 -K + (1 + - 2N X ) cos^, (20) 

where the functions Qq and Qi only depend on the actions I + , K + , and Q 2 is constant : 

Q = -3K+ + 21+ + U 3 l\ + U 2 K 2 + - V 23 l+K + , 

fii = -2U 2 K+ - AU 3 I + - V 2 iK + - V 31 I+ + V 23 (I + + 2K+), 

Q 2 = U x + U 2 + AU 3 + V 21 + 2V 31 - 2V 23 . (21) 

Numerically Q2 is given by ^2 — —3.9642. 

We notice that Hamiltonian (|2*U|) does not depend on the angles 8 2 , and 9 3 which is equivalent 
to saying that I + and K + are invariant (constants of motion). 

In order to visualize the dynamics, we perform numerical integration of the equations of 
motion associated with Hamiltonian (1201) . 

A plot of some trajectories for J + = 1 and different values of I + in the accessible region of 
the phase space (grey region) is given in Fig. El In these plots, we have transformed from 
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action-angle (N,9i) to cartesian coordinates Pi = \fN\ cos 9\, Qi = y/Nisin9i. In these 
last variables, the accessible region is, for K + < 0, i.e. 1+ < JJ+/2, the set of points in 
the interior of the circle Aq = I+/2 and such that these points are inside T\ or outside T\ 
(where the canonical transformation from the original variables bj to action-angle variables 
(Nj, 9j) is well defined). For K + > 0, i.e. I + > J+/2, the accessible region is the annulus 
K + < Ni < I+/2 and such that these points are inside T 2 or outside T 2 . 
We clearly see that these plots are characterized by elliptic, hyperbolic points and the 
associated separatrices. Below, we discuss the location of these points with respect to the 
parameter I + . 

The fixed points are obtained when i\q = and 9\ = 0. As depicted in Fig. ^ there is three 
families of fixed points. The first family, denoted e(I+), is located at 9\ = and is elliptic. 
Its position Aq is solution of the equation: 



W 

fii + 2fi 2 Ax + — 



'Aq - 1K_| 



N 



Ni — K, 



(I+-2N 1 )-4JN 1 (N 1 



0. 



(22) 



which is the equation of 9\ = for 9\ = 0. 

There exist two others families of fixed points located at 6\ = tt, one elliptic point denoted 
ei(I + ) and the other hyperbolic point denoted hi(I + ). 
The location Ni is given by the equation: 



fij + 2VL 2 N X 



W 

Y 



IN - K 4 



N 



+ 



N 



Ni — K H 



(1+ - 2JVi) - AJNiN - K+) 



0. 



(23) 



which is the equation of 9\ = for 9\ — tt. 

The location of all these fixed points as a function of 1+ is depicted on Fig. [T] 
We see on the different plots of Fig. El that when I + is close to zero, there is no fixed points. 
When we increase I + there is only one elliptic fixed point e(I + )( Fig. a) ) which appears 
near zero. As the value of I + increases, there are three fixed points: besides the previous 
elliptic fixed point e(I+), there is a bifurcation (fold, which is structurally stable) where two 
new fixed points appear, an elliptic point ei(I+) and a hyperbolic one /ii(I + )(Fig. Ofb)). 

We have a homoclinic connection for the hyperbolic fixed point. When the value of I + 
still increases and remains lower than J+/2, the elliptic fixed point e(I + ) moves to the right 
and the hyperbolic one fta(I+) moves to the left, both going towards the border of the phase 
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space. At the same time, the elliptic fixed point ei(I+) moves to the right going to the 
center of the phase space so approaching Ni — 0. When I + is close to J+/2, the elliptic fixed 
point ei(I + ) remains near the border iV\ = but it is not apparent in the figures and the 
hyperbolic one hi(I + ) remains near the border Ni = I+/2. (Fig. E(c)). 
For I + > J+/2 the hyperbolic point hi(l + ) remains at the border N% = I+/2 and the elliptic 
point ei(I+) is in the border of the annulus where K + > 0. At the same time the region 
in which N% evolves decreases ( Fig. Efd) and Fig. HHe)). For 1+ near 2J + , the phase space 
merges to the point S. 

We notice that for I + = J+/2, the resonance occupies a larger region of phase space than the 
other cases. This is in agreement with Fig. ^ All these fixed points correspond to periodic 
variations in variables Ni,N 2 ,N 3 and $ 1 ,$ 2 ,$3. Indeed, due to the invariants given by 
Egs-tlSjl- lfrUj) . if N\ is fixed, N 2 and N 3 remain fixed and the angles $i, $2, $3 evolve linearly 
in time. For example in the case I + = J+/2 ( i.e. K + = 0), we give explicitly the dynamics 
corresponding to e(J+/2). The equations of motion are 

= V^i(I + -2AMsin0i, 

at 

Jj± = Q 1 + WI+cose 1 + 2N 1 (Q 2 -2Wcos0 1 ). (24) 

For 6*i = and 

Nl = -2Q 2 + 4W> (25) 
we get an elliptic fixed point with Ni < I+/2 for all I + = J+/2. This means that this fixed 
point exists for all admissible values of the parameter J + and for I + = J+/2 . 
We deduce the expressions of the linear motion of the angles : 



$i(t) = l-V 31 I + + -WI + + (2U 1 + V 21 -2V 31 -W)N 1 t + $i(0), 
= c\t + an, 

$ 2 (0 = 3 + V23I+ + ^WI + + (V21 + 2C/ 2 - 2V23 - W)Ni t + $ 2 (0), 
= 9c 2 t + a 2 , 

$ 3 (t) = [2 + 2U 3 I+ + (-V 31 + V 23 - AU 3 + W)Nx] t + $ 3 (0), 
= 4c 3 t + a 3 . 
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This linear dynamics gives the following shape of the wave: 



r)(x,t) 




N\ cos (x + C\t + ai) 



+v / 3 v // iVi cos [9(x - c 2 t) + a 2 ] 



+ y/2 yjl+ - 2N l cos [4(z - c 3 t) + a 3 ] 



(26) 



which is a sum of three traveling waves with different velocities. 

For this case, we have plotted on Fig. 01 the elevation of the surface rj(x,t). We see that r) 
is quasiperiodic in time. 

D. Exchanges of energy on C+(I+) 

We will now see how the exchanges of energy occur among the three modes for different 
values of I + . 

Since all the orbits of C+(I+) are closed, the function Ni(t) is periodic. From the two 
invariants I + and JJ + , the two other actions N 2 (t) and N 3 (t) are periodic with the same 
period as Ni(t). Therefore, all the exchanges which occur in the system are periodic in 
time. 

• We consider C + (I + ) with 4±- < I + < 2J + . Let us first consider the trajectories with 
initial conditions close to S, i.e., such that iV 2 (0) and N 3 (0) are of order e « 1. If 



In fact, this behavior is valid for all values of I+. So, if we have initially N\(0) >> N 2 (0) 
and A^(0), then by looking at the invariant I + , we see that for A^(0) close to zero, Ni(t) 
decreases and at the same time N 3 (t) increases. The invariant 2N 2 + N 3 w 2N 2 (0) 
shows that N 2 (t) also decreases and N 3 (t) does not grow more than 2A^(0) which 
is small as we can see in Fig. where we have plotted the time evolution of the 



we set I + = 2JI+ - e, we get from Eqs. ©-(HDD 



!±_£ <JVl(i )<!±, 

2 6 ~ W ~ 2 ' 



< N 2 (t) < - 
6 

< N 3 (t) < |, 



for all t. 
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amplitudes Ni, N 2 and A/3 for the initial condition (iVi, #1) = (-±, 7r). 
Therefore, when almost all the energy is initially carried out by the lower frequency 
mode N\, the oscillations of the amplitudes are very small and there is no significant 
energy exchange between the modes. 

By looking directly at Fig.^ we also see that, for such values of I + , only the modes N 2 
and N 3 may exchange their energy while N% may only sustain small energy fluctuations. 
In the same way, for an orbit corresponding to C + (I+) with < 1+ < JJ+/2, there is 
no significant energy exchanges for initial conditions close to P, and in any case only 
N\ and A/3 enable some energy exchange among the modes. 

• In order to obtain large energy exchanges between the three modes, we have to put as 
much energy in the side bands as in the central carrier: the contribution of energy to 
the side bands is equal to ooiNi + uj 2 N 2 = Ni + 3N 2 . The maximum energy exchange 
is obtained when N% + 3N 2 = 2N 3 = JJ+/2. This happens in the zone where N\ and 
N 2 are equally important namely close to C + (J + /2). 

In this case, using the form of the three integrals of the motion and assuming I + = 
J+/2 + e, we get: 

0<N 3 (t)< 3 -± + e, f<^i(*)<Y + |- W 

Then if we assume that at t = AT 3 (0) >> iVi(O) and N 2 (0) such that the most energy 
is stored in the highest frequency mode, from 2N\ + A3 ^ AT 3 (0) and N\(t) — N 2 (t) = 
iVi(O) — iV 2 (0) we see that N 3 (t) decreases and at the same time Ni(t) and A^ 2 (t) 
increase. So nothing prohibits that energy will be completely transferred from the 
highest frequency mode into the low frequency modes. What we observe is a periodic 
modulation of the wave amplitudes. We see that most of the initial energy carried out 
by the central frequency mode of the wave, A/3, is transferred to the side bands and the 
magnitude of the amplitude oscillations is very large. This is the usual Benjamin- Feir 
instability: the carrier and the two side-bands evolve on a periodic cycle of modulation 
and demodulation known as Fermi-Pasta-Ulam oscillations. In Fig. H we see how 
the energy exchange is done among the amplitudes Ni, N 2 and N 3 for for initial 
(N 1 ,9 1 ) = (I + /8,0). 

For this initial condition, we have plotted in Fig. the shape of the wave f]{x,t) for 
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this initial condition. The periodic modulation and demodulation of the three modes 
give a shape of the wave which is a modulated traveling wave. 

• We can also see exchanges of energy for I + = J+/2 i.e. on the segment C + (J + /2) and 
for initial conditions so that A 3 (0) is close to zero and A 1; N 2 lies on C + (J + /2) so 
that JVi(O) = iV 2 (0) « J+/4. We recall that if A x (0) = N 2 (0) then Ni(t) = N 2 (t) and 
K+ = 0. 

For this case, from Eq.{(SJ) , we conclude that Ni(t) and N 2 (t) decrease and N 3 (t) grows 
until it reaches its maximum which is J+/4. As shown in Fig- [HJ the trajectory stays a 
long time close to its initial position and then make a fast excursion before an equally 
fast return to the departure point. The fast oscillation corresponds to an important 
exchange of energy between the carrier and the side bands since the carrier almost 
reaches its maximal energy during the burst. We see in the Fig. EI the shape of the 
wave in this case which is quasiperiodic. 

IV. DYNAMICS WITH TWO RESONANCES 

In this section we analyze the dynamics of Hamiltonian (jHJ with the two resonances. 
The phase space is now the product of E + by E_ associated with each resonance. In this 
case, the dynamics shows a rich variety of regimes allowing a more efficient energy transfer 
among the modes. We give the dynamics of different orbits in the phase space for different 
values of the parameters and also give the type of the shape of the wave r](x,t). 
As in the case of one resonance, we start by simple situations: we explore the dynamics of 
the points P + x P_, S+ x S-, T + x T_ and the segment (P+S+) x (P_S-) for which we get 
the analytical expression of the wave 7](x,t). 

A. Dynamics on points P + x P_, S+ x S- and T + x T_ 

Let us start with the case where the position of the projection of the orbits on (E + ) x (E_) 
is fixed at P + x P_ which correspond to I + = I_ = (see Fig. Since I + = 2N\ + A 3 is a 
conserved quantity and N±, A3 are such that Aq > and A3 > 0, we have Ni = A3 = for 
all time and A 2 = J+/3. Similarly we have A_i = A„ 3 = and A_ 2 = J+/3. 
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Then for P + x P_, all the actions are constant and the dynamics is connected to the evolution 
of the angles. The corresponding dynamical equations are: 

.db 2 (t) 



dt 
. db_ 2 (t) 

dt 



The solution in bj(t) is then : 

b-2(t) 

Jo + I- - 2 J 4 

a 



(3 + 2U 2 (N 2 - 2N_ 2 ))b 2 , 
(3 + 2U 2 {N_ 2 - 2N 2 ))b_ 2 . 



(28) 



(29) 



A /j+/3exp(-z(9ct + 0+)) , 
y{j+/3 exp (— z(9ctct + </>_)) 
0, for j = 1,-1,3,-3, 



where 



1 



'1 



J+ - 2J_ 



Jo + J+ - 2J_ 3 
We notice that a is related to J + and JJ_ by 



and 



(1 + 2a)J_ - (a + 2)J + = (a - 1)JJ„. 



The expression of the shape rj(x, t) is given by: 



r)(x,t) 



-1/4 



—= £ l^'l 174 (% + aI i )^ + (a* + a_ i )e-^ 
2vrV2 i=1)2>3 V 



(30) 



(31) 



(32) 



(33) 



where the variables aj, as function of &j are given by (see Remark 2 of the Appendix): 

>>, ■ 3i};;;; ; . for j = 1,2,3. (34) 

Similarly a_j is given by replacing kj by — kj in Eq. ()34jl. and a, = for j ^ 2, j ^ —2. 
Then we deduce the surface shape: 

-1/4 



7?(x,t) 



9 



7T 



>/2 



+ cos(9(x — ct) — 4> + ) + y J- cos(9(x + act) + </>_) 

+1/1-^1+ cos(9(a; + (2a + l)c£) + 20_ + <p + ) 
+VI + 7jTcos(9(a; - (a + 2)ci) - 2<f) + - <f>_) 



(35) 



where V = V$$_ 2 - 2 « -0.082. 



The wave is the sum of four traveling waves with four a priori different velocities 
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c, ac, {2a + l)c, and (a + 2)c. 

According to the values of J+ and J_ and more precisely, depending on a, these solutions 
are fixed points, traveling or standing waves, periodic or quasiperiodic solutions. 



If JJ + = JJ_, then a = 1. We deduce the shape of the wave : 
9 



r}(x,t) 



- 1 / 4 v / 2 



J + cos 9x 



71 



x ( cos(9ct + 



+ VJ+ cos(27ct + 3- 



>+ + 



2 2 
which is a standing wave. Such wave is represented in Fig. HUf a). 



For the particular case J + = J_ = J , a is not defined and the velocities c and etc 
vanish. The expression of the surface wave is: 



rj(x,t) 



9' 



7T 



X COS 



cos(9a; 



4>+ - 4>- 



+ 



-) + cos(3- 



2 2 

The wave is stationary (see Fig. ITUT b)) and its amplitude is a function of the initial 
phases. 



When J + + J_ = 2 J so that a = —1, we have the unique family of traveling waves 
with velocity c whose shape is given by: 

7T\/2 

+VJ_ v /j7cos(9(x - ct) + 



r](x,t) 



cos(9(a; — ct) — <p + ) + -i/JL cos(9(x — ct) + 0_) 



+ 



+VJ + V /J_ cos(9(x - ct) - 20+ - 0_) 
Fig. ITUT c) shows the elevation of this surface. 



(36) 



Time periodic waves are obtained for a G Q including the particular case — J + + 2 J_ 
Jo and — J_ + 2J + = J , for which a is respectively indefinite or zero. 
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The ratios between the four angles, evolving in time, which are —ct, act, (2a + l)ct 
and — (a + 2)ct are rational. The wave is periodic both in space and time, forming 
the square-shape patterns as observed in Fig. ITUT d). 

• In the case where a is irrational, the wave is quasiperiodic in time, expressed as a 
sum of four traveling waves with different velocities. We have plotted such wave in 
Fig. EHe). 

We notice that a similar variety of regular surface shapes ( standing, stationary, traveling, 
time-periodic and quasiperiodic waves) are obtained for the two other points of (E + ) x (E_) 
which are S + x 5_, T + x T_, and for the points P + x SL, P + x T_, S + x T_, S+ x P_, 
T + x S-, and T + x P_. 

B. Dynamics of a point on the segment (P + S + ) x (P 5* ) 

We consider the point which is at the intersection of C+(I+) x C_(I_) and (P + S'+) x 
(P-S-). This point is such that 63(0) = and 6-3(0) = 0. From equations of motion of 
Hamiltonian (jSJ), we deduce 63 (t) = and 6-3 (i) = for any time. From the invariants 
we obtain that N±, N2, iV_i and A^_2 are constant and equal to Ni(t) = I+/2 and A^(t) = 
(2J + - I + )/6, = I_/2 and AL 2 (t) = (2J_ - I_)/6 . The equations of motion become: 




We deduce the expression of the shape of the wave : 




+yl- cos(x + c_it + a_i) + a/2J_ — I_ cos(9(x — c_ 2 t) + a^ 2 ) , 





(37) 
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where c x = 1 + (U x + V 21 /6)I + - {2U X + K- 2 i/6)I_ + V 2 iJ + /3 + V_ 2 iJ_/3, 
c 2 = [3 + (-U2/3 + ^2i/2)I+ + (2Z7 2 /3 + VC 2 i/2)I_ + 2f/ 2 J+/3 - 4C/ 2 J_/3] /9, 
c_i = 1 - (2^ + F-2i/6)I+ + (f/i + K 2 i/6)I_ + V_ 21 J+/3 + V 2l S-/S, and 
c_ 2 = [3 + {2U 2 /3 + V-. 21 /2)I+ + (-E/2/3 + ^2i/2)L_ - 4C/ 2 J+/3 + 2C/ 2 J_/3] /9. 
This wave is the sum of four traveling waves with four different velocities. 

C. Dynamics of orbits on the segment (P + S+) times the interior of C_(I_) 

We consider the dynamics of (P + S + ) times the interior of C_(I_) (or similarly the 
interior of C+(I+) times (P_S'_)). We have 63 (t) = for all t, so AT 3 (t) = is a constant of 
motion and Ni(t) and N 2 (t) are constant for all t. We have now a system with five degrees 
of freedom bi, b 2 , 6-2 and 6_3 with five constants of motion, I + , I_, JJ + , J_ and H. 
It is then integrable. The modes b\ and b 2 evolve with a constant modulus and a phase 

depending on |6-i|, \b- 2 \ and 1 1> 3 1 . The dynamics of the modes b-i, b- 2 and 6-3 in C_(I_) 

depend on the second resonance just through the modulus of b x and b 2 which are constant. 
So the exchanges of energy on C_(I_) are not influenced by the evolution of b\ and 6 2 . The 
dynamics of the modes b-i, &- 2 and 6_3 is then similar to the one studied in the case with 
one resonance in Sec. IIHI 

D. Dynamics in the interior of C + (l + ) x C_(I_) 

In order to study the dynamics on the nontrivial sets C + (I + ) xC_(L) on the interior of the 
phase space (E + ) x (E_), i.e. such that iVj- are strictly positive for j = —1, —2, —3, 1,2,3, 
we perform a canonical transformation in action-angle coordinates, bj = y'/Vje -1 ^ and 
Hamiltonian (jHJ) becomes: 

H= E ^•(iV i + AL i ) + E Uj{Nf + Nlj - ANjN-j) 

.7=1,2,3 i=l,2,3 

+K 21 (AT 1 AT 2 + AT_!iV_ 2 ) + VC 2 i(iViiV_2 + N_iN 2 ) 
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+V 2 3(N 2 N 3 + AL 2 AL 3 ) - V- 23 (N 2 N- 3 + N- 2 N 3 ) 
+W^N 1 N 2 N 3 cos($i + $ 2 - 2$ 3 ) 

+W^nZn7 2 N^ 3 cos($_i + $_ 2 - 2$_ 3 )- (38) 

Using the invariants I+, I_, J + and J_ , we now proceed to the reduction of the dimension 
of the system to two degrees of freedom. 

We perform as in the case with one resonance the canonical transformation from 

(N u N 2 , N 3 , AL 1; N_ 2 , N_ 3 , $ 1; $ 2 , $ 3 , $_ 2 , $_ 3 ) to 
(iV\, X+, I+, i\Li, X_, I_, 0i, 3 , 0_i, 0_2, 0-s): 



iVi 


= N u 


K + 


= N, - N 2 , 


1+ 


= 2iVi + N 3 , 


0i 


= $i + $ 2 - 2$ 3 , 


02 


= -$2, 


03 


= $3, 


iV_i 


= AU 


K_ 


= iV_i - AL 2 , 


I_ 


= 2iV_! + iV_ 3 , 


0-i 


= $_i + $_ 2 - 2$ 


0-2 


= -$-2, 


0-3 


= $-3- 



The actions Ni, K+, I+, ALi, K_ and I_ are now respectively conjugate to angles 0i, 02, 
03, 0_i, 0_2, and 0^3. In what follows, we use N±, iV_i instead of Ni and iV_i for simplicity. 
Hamiltonian (J3%|) becomes: 



H = tt + fliJVi + fi_iiV_i + Q 2 {N% + Nlj) + tt^N. 
+wJn[Jn 1 -K+(I+ - 2iVi) cos 0i 



+WJN_ 1 JN_ 1 - K_(I_ - 2ALi) cos0_i, 



(40) 
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where the functions Qq, f2i and Q-± only depend on the actions I+, K+, I_, K_ and f2 2 and 
f2 3 are constant : 

Q = -3(K_ + K+) + 2(I_ + 1+) + U 3 (f_ +I 2 + - 4I + I_) 

+C/ 2 (K?. + K + -4K + K_) 

-V r 23 (I+K + + I_K_) + V1 23 (I+K_ + I_K+), 
fij = 2C/ 2 (-K + + 2K_) + 4C/ 3 (-I + + 2I_) - V 2 iK + - VC 21 K_ 

-V 31 I + + ^-3il- + ^ 23 (I+ + 2K+) - V_ 23 (I_ + 2K_), 
= 2C/ 2 (-K_ + 2K+) + 4C/ 3 (-I_ + 21+) - V 2i K_ - ^_ 2 iK+ 

-y 31 L + V_ 31 I+ + V2 3 (L + 2K_) - V_ 23 (I+ + 2K+), 
fi 2 = U x + U 2 + 4C/ 3 + V21 + 2V 31 - 2K 23 , 

fi 3 = -AU X - AU 2 - I6E/3 + 2V_ 21 - 4V1 31 + 4V_ 23 . (41) 
Numerically f2 2 and f2 3 are given by: f2 2 = —3.9642, fi 3 = —16.0933. 

Since Hamiltonian ()40|) does not depend on # 2 , 3 , 0_ 2 and 0_ 3 , we can consider I + , K + , I_, 
K_ as parameters of the system for the dynamics of Nx, iV_i, 0i and 

We recall that Ni and iV_i are such that max(0, K+) < Ni < I+/2 and max(0, K_) < N_i < 
I_/2. The accessible part of phase space is thus contained into the product of two interiors 
of a circle or an annulus, depending on the values of I + and I_. 

In order to visualize the dynamics on C+(I+) x C_(I_) we compute Poincare sections of the 
dynamics and then we project these surfaces on (Ni, 0i)-plane. The Poincare sections S are 
constructed in the following way: 

£ = I (iVi, iV-i, 01, ^-i) : 0-i = 0, < 0, # = (JV ls JV_ ls ^ , 1 + , . (42) 

Now the parameters are the energy H, and the four invariants I + , I_, J + and J_. 

Since Hamiltonian (|3*Hj) has two degrees of freedom, Poincare sections are two dimensional 

and therefore we get an accurate picture of the dynamics by looking at the projections of S 

on (N\, 61). For most of the values of the parameters, these sections show a rich variety for 

the dynamics, a mixed phase space with regular and chaotic behaviors. 

We give a plot of some projections for J + = J_ = 1, H — 2 and for different values of 

1+ = I_: from !+ = !_ = 0.3 to 1+ = I_ = 0.8. 

We observe on Fig. [TIT a), when I + and I_ are near zero so that we are in the neighborhood 
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of the point P + x P_, the motion is completely regular and all the orbits lie on invariant tori. 
For this value of H, the dynamics occurs close to the center of the accessible region of phase 
space. The elliptic and hyperbolic fixed points of the Poincare sections are of the same type 
as those of the elliptic and hyperbolic orbits of the case with one resonance. Besides these 
points, there are one elliptic and two hyperbolic points due to the second resonance. We 
recall that the invariants depend on the variables of only one resonance. So, for small values 
of I + and I_ we can make the same reasoning as for the one resonance case and we can 
conclude that in the neighborhood of the point P + x P_, the exchanges of energy between 
the modes are not significant and it leads to a regular behavior of the system. Thus in the 
system with small values of I + and I_, the only possible motion for the system is periodic 
or quasi-periodic in time, corresponding to regular shapes for the wave r)(x,t) (see Sec. 1111)1 . 
To undergo a transition to chaos, we increase the parameter I + ( which increases the effect of 
the nonlinear terms up to a certain threshold), and for a value of I + in the neighborhood of 
I + = J+/2, we see in Fig.lTlTb) and Fig. lTlT c) , that the phase space is made up of regular and 
chaotic motions. Chaotic trajectories are located between regular trajectories (on invariant 
tori). We observe that the motion is clearly affected by the vicinity of hyperbolic points: 
the irregular motion appears in the neighborhood of hyperbolic periodic orbits breaking up 
their neighboring KAM tori. 

The first of the two regular zones which surround the chaotic sea (outer region), is obtained 
for large amplitudes of the modes of low frequency Aq and A_i which are larger than the ones 
of high frequency modes A3 and AL3. The second regular zone (inner region) is obtained 
when the amplitude of the modes of low frequency Aq and A_i are small. In these regular 
zones, the exchanges of energy are not significant and give rise to regular behaviors of the 
system. On the other hand, for orbits close to the hyperbolic points, the modes oscillate 
chaotically (see Fig. IT!?)) . As the parameter I + increases, the region in which Aq evolves 
shrinks and all the orbits become stable, the system admits only regular behaviors as it can 
be seen on Fig. fTTT d) . This is due to the fact that in the accessible region, the values of 
Aq are larger than those of A3 and for this case the exchanges among the modes are not 
significant. 

We conclude that chaotic exchanges of energy may be done in the neighborhood of C + (^-) x 
C-(*f). In Fig. [U we give the plot of the modes Nj(t),j = 1,2,3,-1,-2,-3 and in 
Fig- El we give the shape of the wave rj(x, t) for I + = I_ = J+/2 and for an initial condition 
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(Ni,9i) = (0.152,-3.07) which is in the chaotic region. In these two plots, we see the 
irregular motion of the modes Nj (chaotic exchanges of energy) , and the chaotic structure 
of 7](x,t). We notice that in the chaotic zones, significant energy exchanges occur among 
the modes. 



V. CONCLUSION 

Following the framework of Zakharov in Ref. and Krazitskii in Ref. Q], we have built 
a Hamiltonian model for six interacting waves for gravity water waves in infinite depth. 
The interaction are those of two coupled Benjamin- Feir resonances. The knowledge of four 
constants of motion allowed us to reduce the system to a Hamiltonian system with two 
degrees of freedom. We have first studied the dynamical properties of the Hamiltonian with 
only one resonance which is integrable. This simple case allowed us to understand the basics 
of the full dynamics. A whole variety of solutions have been found in this simple case such 
as traveling, time periodic and quasi-periodic waves. We showed that exchanges of energy 
occur in the zone of phase space where the energies initially put into the two side bands 

are equally important and when the energy in the central mode is the same as the sum of 
the energies of the two side bands. This situation corresponds to phase space regions close 
to the separatrices. All these exchanges are done in a periodic way. Thus, when we take 
into account only one resonance, the only possible motion is a regular time quasiperiodic 
(and space periodic) wave. 

Then we have studied the dynamics with the two resonances. Besides the regular waves 
such as traveling, time periodic and quasiperiodic waves, the full system exhibits stationary 
and standing waves. We have visualized the dynamics using numerical Poincare sections. 
We have showed that the system with the two coupled resonances undergoes a transition to 
chaos. Chaotic behaviors occur in the regions where the energies of the side bands of each 
resonance are equally important and near the broken separatrices. In this case, there is a 
strong exchange of energy between all the modes and it is done in chaotic way. 
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DERIVATION OF THE HAMILTON! AN GIVEN BY EQ. © 



In this appendix, we briefly derive the expression of Hamiltonian (JHJ) from the Hamil- 
tonian derived by Krazitskii in Ref. [y| given by Eq. (2.11). We get the expression for 
the Hamiltonian in the case where the three modes (denoted by 1, 2 and 3) are such that 
k\ + &2 = 2/^3 and uo\ + c<j 2 = 2u^. 

In this case, we notice that there are no terms of order 3, i.e. terms of the form a^a^a k n 
such that k = k' + k" for k, k\ k" G {k\, /c 2 , A; 3 , —k\, —k 2 , —ks}. 

We expand this Hamiltonian up to order four, i.e. we neglect terms of order five in Hamil- 
tonian (2.11) of Ref. since we are in the case of weak nonlinearity ( small amplitudes of 
the wave). Hamiltonian (2.11) of Ref. Q] becomes: 

H = U i( a i a *i + a ~i a *-i) + 

i=l,2,3 

+6 ^ V^_jAa*_ id-ia-jdj + a*aiaja_j + c.c) 

i,j=l,2,3, 

+3V2 _i 33(Q2 a -l a 3 + CL*-2 a l a -3 + C.c) + 3^/,-2,3,3( a l a -2 a 3 + 1^2^-3 + C.c) 

+SVsZs ! i,2( a 3 a -3 a i a 2 + al 3 a3a_ia_2 + c.c) 

+9 12 v ffiA a f a2 i + a -i a -i) + 2 J2 v l 2 \i,- l ( a * a U a i a 'i) 

/ 8=1,2,3 i=l,2,3 
+2 ^ ( Vij,i,j( a i a *j a i a j + 0-li a lj a -jO-i) + ^-j.i.-jC * - j a i a -j + a *-i a j a -i a j) ) 



i,j=l,2,3 



+2 V -u,-j,j( a -i a i a -i a i + C ' C ) + ^2,3,3 K°2 a 3 + a-lfl-2 a -3 + c - c ) 

i,j=l,2,3, 

+2V 3 (2) li _3 j2 (a3a!_ 1 a_ 3 a2 + al 3 a*a 3 a_ 2 + c.c) 

Z 1=1,2,3 »>J=1,2,3, 

+3V r 1 (4 2 ) i _3,- 3 («-l«-2«3 2 + «l«2«-3 + C-C), (43) 

where the coefficients verify the symmetry V^f kl = VZi_ j _ k i } for all n = 1, 2, 4 and 



Q 



fc, Z G {1, 2, 3, -1, -2, -3} as well as the other symmetries mentioned in Ref. |lfl|. We 
verify in the expressions of the coefficients given in Ref. Q] that all the coefficients V i ^_ ii 
are equal to zero. 

The canonical transformation from the variables a k to the new canonical variables bk, 
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that maps Hamiltonian (|43|) to its normal form up to terms of order 5 in bj is given 
at = h + 6^2,2(61626-2 + 6 1 6*6*_ 2 ) + 6Wji-3, 3 (W-3 + WblbU) 

+6^-1,1(6-2^2^! + 6 2 6*6*_ 1 ) + 6V$,_ 1 , 1 (&- 3 &i 3 &-l + 636561i) 
+3y 2 ( ) 1 i,3,3&-2&! 2 3 + 3^2,3,36-2^ + 3^3,1,2636-36; 

+2V^ 2 i!i,-2, 2 6-i6-26 2 + 2Vi 2 i ; i) _3 i3 r i 6- 3 63 + ^V%^ 2 hb*_^ 2 
+3Vl%_ lt _ 1 b%b\ + ^ 4 i,_ 2 , 2 6* x 6* 2 b* 2 + ^,1,-3,36-16-36,; 

+3V^ 1 ^_2,3 ! 3&2^ > — 2 3' 



a 2 = b 2 + 6V^ ) _ 2i2 {b. 1 b*_ 1 b*- 2 + b x b\b*_ 2 ) + 6^-3,3(62636-3 + 62656*3) 
+6^_ 1 ,i(6 2 6 1 6_ 1 + 626^) + 6V > 3 ( i 1 3 ;_ 2i2 (6_36*_36*_ 2 + b 3 b*b*_ 2 ) 
+3^2 ( , 1) i,3, 3 6-i6 2 + 3t/ 1 ( , 1 i,3,36-i6! 2 3 + 3^3,^636* 3 
+2V r i 2 2 !2,-i,i6- 2 6-i6i + 2^2,-3,36-26-363 + 2Vgl h _ 3t2 b 3 b*_ 3 b-i 
+3V 2 ( ,J- 2 ,- 2 6! 2 2 6; + V%^ 2 b*_ x b\b*_ 2 + ^2,-3,36-26-36; 
+3V; (4 2!-3,-36i6! 2 3, 

a 3 = 63 + 6V 1 ( i i;_3 i3 (6-i6*_ 1 6*_ 3 + 6^6*3) + 6^ 2 ( i 1 2 !_3 i3 (6_ 2 6*_ 2 6*_3 + b 2 b*b*_ 3 ) 
+6^S-i,i(6 3 6i6_ 1 + 6 3 6?6* J + 6^2,2(63626-2 + 6 3 6;6*_ 2 ) 
+6W,- ) i,3,36 2 6*-i6; + 6^ 1) 2 ,3,36i6: 2 6; + 3^3^2(6-36261 + 6- 3 &*_ 1 &*_ 2 ) 
+2^ 3 ,-i,i6: 3 6-i6i + 2^3^2,26-36-262 + 2y 3 ( i 2) 1> _3 i 2(6- 3 6*_26i + 6_ 3 6 2 6! 
+3^_3,- 3 6- 2 3 6; + V%^ z b\b\b% + V^^UWU 
+6^?- 2) 3 1 36* 1 6* 2 6*. 
Similarly a_j is given by replacing kj by — kj in Eqs. (j44|) . (|45j) .and (|46|) . 



The coefficients vff^, vff^ m and v}j\ m are given by : 

V. (1) 

t/(1) _ i,j,l,m 



,i,L,m ) 
UJi - UJj -uji- UJ m 

v {2) 

T/(2) v i,j,l,m r 1 / 1 

^' l ' m UJi + UJj -UJi-UJ m 
T/(4) 

y(4) _ "t,j,i,m 
i.j.l.ra 1 1 1 
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Remark 1 : We notice that when the modes of one resonance are equal to zero, e.g, 
6_j = for j — 1, 2, 3 then we have aj = bj for the second resonance (j = 1, 2, 3). 

Remark 2 : When two modes of each resonance are equal to zero e.g bj = b^j = for 
j = 1,3 then a 2 = b 2 + 2V 2 (4 2 \_ 2t _ 2 b%b* 2 , where V 2 { %_ 2 _ 2 = -0.08188. 
The transformation which maps {aj} into {bj} is canonical i.e. it satisfies 

idbj A db* = ^ idaj A da*. (48) 
j j 

This transformation maps Hamiltonian (J43|) into Hamiltonian (jHJ) with the following val- 
ues of the coefficients W and Vij. 







u 2 


u 3 


W 






0.0063 


4.61645 


0.40528 


0.52648 
















v 21 


V-2X 






v 23 


V- 23 


0.07599 


0.53193 


0.02533 


0.17312 


5.16738 


2.12774 
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FIG. 1: Representation of the phase space E + and the three families of fixed points found on 
(Ni,8\) plane. The location of the elliptic points e(I+), ei(I+) is represented by a dashed line. 
The location of the hyperbolic points h\(I+) is represented by a dotted line. The bifurcation 
creating the pair of elliptic/hyperbolic ei(I+)//ii(I+) point is indicated by a dot. 

FIG. 2: Elevation of the surface rj(x,t) corresponding to P (1+ = 0) for Hamiltonian (|12|) with 
one resonance. 

FIG. 3: Phase space trajectories of Hamiltonian (|12|) with one resonance for different values of 
K+: (a) 1+ = 0.1, (6) 1+ = 0.12, (c) 1+ = 0.5, (d) 1+ = 0.505, (e) 1+ = 0.6, (/) 1+ = 1.7., The grey 
region is the accessible region of phase space: which is the annulus max(0,lK + ) < Ni < I+/2. The 
bold curve T\ represents the trajectory which passes by N\ = 

FIG. 4: Elevation of the surface rj(x,t) for Hamiltonian ()12l) with one resonance for J+ = 1, 
I + = J+/2, and for the fixed point e(JT+/2). 

FIG. 5: Temporal evolution of modes iVj for i = 1, 2, 3, for Hamiltonian (|18|1 with one resonance, 
1+ = 1.2, JT + = 1 and an initial condition (N\,6i) = (I + /2,7r). 

FIG. 6: Temporal evolution of modes Ni for i = 1, 2, 3, for Hamiltonian with one resonance , 
J+ = 1, 1+ in the neighborhood of 1+ = JT+/2, and an initial condition (iV"i, ^i) = (J+/8, 0). 

FIG. 7: Elevation of the surface rj(x,t) Hamiltonian 1)18|) with one resonance for JJ + = 1, in the 
neighborhood of 1+ = J+/2, and and an initial condition (Ni, 9\) = (J+/8, 0) 

FIG. 8: Temporal evolution of modes Ni for i = 1, 2, 3, for Hamiltonian 1)18(1 with one resonance , 
J+ = 1, 1+ = J+/2, and an initial condition on the separatrix (iVi, 9\) = (JT+/4, 2.8). 

FIG. 9: Elevation of the surface rj(x,t) for Hamiltonian (|18|) with one resonance for J + = 1, 
1+ = JJ+/2, and for the fixed point (N u 9i) = (J+/4, 2.8) 
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FIG. 10: Plots of r](x, t) for Hamiltonian J5J) with two resonances and for different values of a: (a) 
standing wave a = 1, (b) stationary wave a indefinite c = 0, (c) traveling wave a = —1, (d) time 
periodic wave a E Q (— J + + 2J_ = Jo), (e) time quasiperiodic wave a = \/0.4. 

FIG. 11: Poincare sections for Hamiltonian (|38|) with two resonances for = 2, J + = J_ = 1 and 
for different values of 1+ : (a)I+ = L = 0.3, (6)1+ = L_ = 0.4, (c)I+ = I_ = 0.5, (d)I+ = L = 0.8. 
The grey region is the accessible region where max(0,lK + ) < Ni < I+/2. 

FIG. 12: Time evolution of Nj for j = —1, —2, —3, 1, 2, 3, for Hamiltonian (j38ft with two resonances 
and for I + = L = 0.5 and for an initial condition in the chaotic region (Ni,8i) = (0.152, —3.07). 

FIG. 13: Shape of the wave r](x, t) for I + = L = 0.5 for Hamiltonian (|38|) with two resonances 
and for an initial condition in the chaotic region (N\, 9\) = (0.152, —3.07). 
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